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SUMMARY 


The quasi vortex-lattice method is reviewed and applied to the evalua- 
tion of backwash, with applications in ground effect analysis. It is also 
extended to unsteady aerodynamics, with particular interest in the calcula- 
tion of unsteady leading-edge suction. Some applications in ornithopter 
aerodynamics are given. 


INTRODUCTION 


The quasi vortex-lattice method (Quasi VLM) has been shown to produce 
good accuracy in lifting-surface problems not only for non-f lapped but also 
flapped configurations (ref. 1). In these applications, the only induced 
velocity to be evaluated is the downwash. However, in some other applica- 
tions such as ground effect analysis and wing-jet interaction (ref. 2), it 
is necessary to compute the u-induced velocity (i.e. backwash) in the 
flow field away from the wing plane. It is the purpose of this paper to 
assess the accuracy of such backwash computation. 

One important feature of the Quasi VLM is the accurate prediction of 
the leading-edge suction without resorting to Kutta-Joukowsky relation. It 
is this feature that makes it possible to extend the method to the predic- 
tion of unsteady leading-edge suction in unsteady aerodynamics. This exten- 
sion is also presented below. 
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SYMBOLS 


aspect ratio 

chord length, m(ft), taken as unity. 

leading-edge singularity parameter, defined in eq. (3c) 
induced drag coefficient 


sectional lift coefficient 


amplitude of sectional lift coefficient in unsteady flow 


total lift coefficient 


lift curve slope, per radian 


pitching moment coefficient 


pitching moment curve slope, per radian 


pressure coefficient difference 

time-averaged leading-edge thrust coefficient 

nondimensional height measured from the wing plane, referred to 
the chord length. See fig. 5. 

flapping amplitude 

flapping amplitude at the wing tip 

reduced frequency, defined as wc/2V 
number of integration points, or Mach number 
number of chordwise vortices 


number of spanwlse vortex strips 


leading-edge suction, N(lb), or wing area, m 2 (ft 2 ) 
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amplitude of unsteady leading-edge suction 

t 

I 

nondimens ional backwash, positive downstream 
amplitude of unsteady backwash on the wing plane 
freestream velocity, m/sec (ft/sec) 
nondimensional downwash, positive downward 

rectangular coordinates, with x positive downstream, y positive 
spanwise to the right and z positive upward 

camber slope 

angle of attack, deg. 
amplitude of pitching oscillation 

propulsive efficiency, percent 
nondimensional spanwise coordinate 
nondimensional vortex density 

oscillation frequency, rad/sec 

density, kg/m 3 (slugs/ft 3 ) 

freestream dynamic pressure, N/m^ (lb/ft^) 

■1/=T 

amplitude of flap angle 

sectional lift coefficient in free air 


Subscript : 
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BASIC ANALYSIS IN QUASI VLM 


For simplicity in presentation, consider the two-dimensional (2-D) downwash 
equation: 


w(x) - 4- f Y(x')d 

^ J 0 x-x’ 


( 1 ) 


The integral is first transformed to a 6-integration and then reduced to a 
finite sum through the midpoint trapezoidal rule. It is obtained that 
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Note that eq. (2) differs from the conventional VLM in that sin0 in the for- 
mulation will eliminate the square-root singularities at the edges and the 
vortex densities are directly predicted, instead of the strengths of discrete 
vortices. Furthermore, the control and vortex points are defined by the so- 
called "semi-circle method." This is illustrated in fig. 1. By solving 
eq. (2) with i^O, N Y k ~values can be obtained. Then the leading-edge singu- 
larity parameter C can be computed by taking i=0, i.e., by taking control 
point at the leading edge. The leading-edge suction is then 


irp- 


(4) 
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In three-dimensional cases, the above concept is also applicable by treat- 
ing each chordwise vortex integral in a similar manner as in eq. (2). In this 
case, not only chordwise control and vortex locations are defined by the semi- 
circle method, but also the spanwise vortex strips and control locations. See 
fig. 1. The detail is referred to in ref. 1. The rate of convergence of this 
method is indicated in figs. (2) and (3) for a 45°-sweep wing of AR-2 and 
constant chord. It is seen that the method converges reasonably fast. 


BACKWASH EVALUATION AND GROUND EFFECT ANALYSIS 


In ground effect analysis, it is known that the image vortex system may 
produce strong backwash to decrease the air velocity on the wing great enough 
to be significant. In fact, as a result the wing lift in ground effect may 
be less than the free air value. This unfavorable effect of backwash is par- 
ticularly important in powered-lift aerodynamics and for wings under high lift 
conditions (ref. 3). Therefore, any formulation of ground-effect problem 
without backwash computation is applicable only to small loading conditions 
as analyzed in ref. 4. 

To see the accuracy of backwash computation by VLM and Quasi VLM, con- 
sider the 2-D expressions for the backwash and downwash: 


u(x,z) 


v(x')dx’ 

0 (x-x') 2 +z 2 


w(x,z) = /' 

2 ’ J 0 (x-x’) 2 + z 2 


If y( x ') = V(1 -x' )/x' , the integration can be performed exactly. The results 
are shown in Appendix A. The downwash expression is included here for later 
comparison. In all computations shown below, equal-spacing elements are used 
for the VLM. For the Quasi VLM, again the 9-transformation is applied first 
before using the midpoint trapezoidal rule. Fig. 4 shows that the backwash 
along the chord evaluated by the VLM at both control and vortex points tend 
to be too high, in particular near the leading edge. On the other hand, the 
Quasi VLM gives quite accurate results everywhere at these not too small z- 
values. With these results in mind, both methods are now applied to the fol- 
lowing 2-D linear ground-effect equations: 


- f 1 

2tt K 


r 1 (x-x') -2h.cn 

f y(x') — dx' 

0 (x-x') 2 +4h^ 2 
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y(x) (l+u(x))dx 


( 8 ) 


c t ■ 2 


r 


The corresponding geometry is given in fig. 5. The second integral in eq. (7) 
represents the normalwash due to the image vortex system. The "2ha" term is 
the backwash contribution. The results of computation are shown in fig. 6. 

Two points are of particular interest. Firstly, the linear vortex theory will 
give better results if the mean surface (as used in the mean surface approxi- 
mation of the linear airfoil theory) is taken through the 3/4 chord point. 
Secondly, both the Quasi VLM and VLM predict approximately the same lift, de- 
spite the fact that the VLM gives higher backwash. This is because the VLM 
also produces higher downwash at a given z-value. This is shown in fig. 7. 

This means that the higher upwash from the image vortex system as predicted 
by the VLM tends to compensate the effect due to the higher predicted back- 
wash. With flap deflected, the VLM predicts lower lift than the Quasi VLM 
does as shown in fig. 8. 

The analysis with eqs. (7) —(8) becomes increasingly inaccur ate at hf0.2. 

As an extreme example, let z=0.05 in eq. (5). With y- V(l-x) /x , the results 
are shown in figs. 9 and 10. It is seen that the backwash is underpredicted 
at the control points and overpredicted at the vortex points by both the VLM 
and the Quasi VLM. Depending on a, this would result in small or even nega- 
tive c^ in eq. (8). The backwash computation with small z has important ap- 
plications in wing-jet interaction theory (ref. 2). Therefore, it is desirable 
to find a practical way to improve the accuracy of the computation. Even 
though increasing the number of vortex points (i.e. the number of integration 
points) can increase the accuracy, it is not a practical way because the number 
of unknowns to be solved would greatly increase, in particular, in 3-D appli- 
cations. Since the inaccuracy is mainly due to the second-order singularity 
in eq. (5) as z-K), a practical method is described in ref. 6 to deal with this 
by weakening the singularity. According to this method, eq. (5) is evaluated 
as follows: 



z 
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where x may be the control points (eq. 3a) or the vortex locations (eq. 3b) and 
Xj - \ (l-cos0j ) , flj - ~ . (11) 
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with M chosen to be 2% for lnterdigitatlon between control and vortex points, 
p being any integer. Note that the last integral in eq. (9) can be integrated 
exactly in 2-D case. However, the similar integral in 3-D would b^ too compli- 
cated to integrate. The accuracy of eq. (10) is also demonstrated with p-3 in 
figs. 9-10. Its application in a wing-jet interaction theory is discussed 
in ref. 7. 


PREDICTION OF UNSTEADY LEADING-EDGE SUCTION 


As mentioned earlier, the Quasi VLM predicts the leading-edge suction 
through the computation of the leading-edge singularity parameter. This fea- 
ture can be easily extended to the unsteady aerodynamics if the downwash ex- 
pression similar to that due to a steady horseshoe vortex can be derived for 
unsteady flow. This has been done recently by integrating the doublet poten- 
tial by parts. The results are given in Appendix B for planar configurations. 
Note that if the oscillating frequency is zero, the expression is reduced to 
that for a steady horseshoe vortex. Using this expression, it is possible to 
extend the steady Quasi VLM to the unsteady case. In 2-D flow, this has been 
done in ref. 8. Some comparison with exact solutions of unsteady leading-edge 
suction (ref. 9) is made in Table I. It is seen that the accuracy of the Quasi 
VLM is quite good. Other aerodynamic characteristics at low or high subsonic 
Mach numbers can also be predicted accurately, including gust response (ref. 8). 

In the 3-D method, note that the downwash expression given in Appendix B 
involves two types of integrals. The integration associated with I-integral, 
eqs. (B.15) and (B.18), can be performed by approximating the integrand, as 
has been done in ref. 10. On the other hand, the integrals, and F^, can 

most conveniently be integrated by approximating the integrands by quadratic 
functions of the integration variable as has been done in ref. 10. The re- 
maining aspects of the method follow closely the steady version. 

To show the 3-D applications, the characteristics of a rectangular wing 
of aspect ratio 2 undergoing the first bending mode of oscillation (ref. 11) 
are computed. To indicate the rate of convergence of the method, the pre- 
dicted complex lift coefficient is plotted in fig. 11 against the number of 
spanwise strips. It is seen that the method converges relatively fast for 

this planform. For instance, with N “4, C is changed by only 1% as N is 

c L s 

increased from 10 to 20. Furthermore, the effect of N is seen to be small 

c 

for this wing without chordwise deformation. The predicted pressure distribu- 
tions at one spanwise station are compared with experimental data and those 
predicted by the Doublet-Lattice Method (DLM) (ref. 10) in fig. 12 with good 
agreement. The 3-D exact solution of unsteady leading-edge suction is not 
available. However, Bennett (ref. 12) has applied Reissner’s high aspect- 
ratio theory to the computation of propulsive efficiency of ornithopters. 
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Here, the propulsive efficiency is defined as 


q SCLV 

— *oo 1* 

n - — 

<l„SC L f, 

__ 1 

c^h - h(n)c £ (n)dn 

-l 

All quantities in eq. (12) have been averaged over one cycle. The i ults for 
two rectangular wings of aspect ratios of 6 and 12 performing line. •• apping 

(h(n) * h t n) are shown in fig. 13. It appears that the high aspect . * o 

theory tends to predict higher thrust than the present method. The agreement 
of the predicted efficiency by both methods is good for AR»12. But for AR*6, 
the high aspect-ratio theory predicts lower efficiency, presumably because it 
would predict much higher input power which depends on the wing loading. 


(12) 


(13) 


CONCLUDING REMARKS 


The quasi vortex- lattice method was shown to possess good convergence 
characteristics in steady wing theory. Its application to the computation 
of downwash and backwash away from the wing plane in 2-D flow showed better 
accuracy than the conventional VLM with z not too small. When z is small, 
both methods become inaccurate in backwash computation. An improved method 
for the Quasi VLM was presented. The Quasi VLM was also extended to the un- 
steady aerodynamics, with the calculation of unsteady leading-edge suction 
being of particular interest. 
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APPENDIX A 


Exact Integration for Eqs. (5) and (6) 


(1) Eq. (5) 



u(x,z) * 2 * J V X' 
0 

+ ~ sign[bx+a 2 ]| 


dx* 
(x 


^ Si J » , lgn[ |, (1 .x ) .,2| 

:-x') 2 +z 2 2,1 *VbV 1 2 


where 


b - |(l-2x) 


i 2 « >/[z 2 +x(l-x) ] 2 +(l-2x) 2 z 2 - [z 2 +x(l-x>] 

2 


(2) Eq. (6) 


w(x,z) 


i~ f 

2 IT J 


(x-x') 


dx' 


(D-D -1 ) cose 2 +D 2 - COS20J 


0 (x-x') 2 +z 2 


D 2 +D _2 -2cos29, 


(A. I) 


(A ?) 


(A. 3) 


where 


n«V2 /{ [23-Z'+(Z' 2 +16z 2 ) is ] lj + [Z , +(Z' 2 +16z 2 ) ii ] i5 | 

1 I 1 — ; ~0 TT^cose) 

l 2+Z + (Z'2+16 z 2 )^ / 


6^ * cos 


Z' * 4z 2 -sin 2 B 


cos0 “ l-2x 


(A. 4) 
(A.S) 

(A. o) 
(A. 7) 
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APPENDIX B 


Downwash due to an Oscillating Horseshoe Vortex 


The downwash produced by an ’’oscillating horseshoe vortex” for a 
planar configuration can be written as 


ii . !!i + 

3z 3z 3z 
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3z r l 2 


— - = F + F 
3z 3 4 


AC f i i 
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2 (y-y ’} 2 


x ».„ (x'-x)(x 2 -x )+e 2 (y 7 -y.)(y’-y) 3 
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^((x 2 -x 1 )(n-y)-HJ)(x 2 -x 1 ) + S 2 (y 2 -y 1 ) 2 (n-y) 


y ,M(y-n) . ,M X 0 M 2 * X 2~*l| r u MR-M 2 x 
« $2 R 8 2 y 2‘ y l V 6 2 


1 ^ 2 “ y 1 )dT (B-5) 


R = v/at 2 +Bt+C 


A = (x 2 -x x ) 2 + B 2 (y 2 -y 1 ) 2 
B * -2[(x-x 1 )(x 2 -x 1 ) + 8 2 (y-y 1 )(y 2 -y 1 )] 


C = (x-Xj^) 2 + e 2 (y-y 1 ) 2 
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(i- 7 = === ) e ' 1 v (T ’% ) dT’ 
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-x 0 +MR 
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(B.16) 


fj 2 = (y-n) 2 


(B. 17) 


^ T l r lH = T l r l at T=0 


( T i r l ) 2 * T 1 r 1 at T= 1 


r^ 2 = rj 2 at t=0 


r^ 2 2 = r^ 2 at t=1 


: 01 = x 0 at T=0 


x n „ = x at t= 1 
02 o 
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+ (1 - — ~ , )e V ( 1 lV( — - - — )- -- - 1 
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B 2 R g 2 y 2 -y l 


(y-Ti) 


x » _i^( T ' +x ) 

e V 0 ; dt' 


T l r l 


(x' 2 +r^ 2 ) 3 / 2 


(B. 18) 


In the above expressions, (x^, y^) and (x^, y ^ represent inboard and 
outboard endpoints of the "bounded element," respectively. 
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Table I. — Comparison of Predicted Leading-Edge Suction 
Parameter with Exact Solutions in 2-D Incom- 
pressible Flow. 

s ■ lim u(x) J 1+x 
X"*'— 1 


s — = 3s/3a , S-? * 3s/36 
a o 


(A) Pitching motion about midchord 


53 — 


k 

N 

Present Method a 

Exact, Ref. 9 

0.1 

10 

1.18870-0.255551 

1.18870-0.255561 

0.5 

15 

0.89889-0.355271 

0.89889-0. 35529i 

1.0 

15 

0. 83374-0. 46737i 

0.83378-0.467481 


(B) Flap Rotation with flap-chord ratio * 0.3 


k 

(N,, 

N,) 

6 

Present Method 

Exact, Ref. 9 


X 




0.1 

(12, 

7) 

0.37115-0.146961 

0.37139-0.147051 

0.5 

(12, 

7) 

0.17483-0.115651 

0. 17496-0. 11576i 

1.0 

(12, 

7) 

0.12972-0.065801 

0. 12982-0. 06592i 

Note: 

N l = 

number 

of doublet elements on 

the airfoil. 


n 2 « 

number 

of doublet elements on 

the flap. 
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SPAN WISE DISTRIBUTION 

Figure 1.- Illustration of chordwise and spanwise control 

and vortex locations. 



Figure 2.- Effect of vortex-lattice arrangement on aerodynamic 
characteristics of a 45°-sweep wing of AR * 2. 6^ ■ 0°; 

M * 0. 
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Figure 6.- Comparison of computed airfoi 1 lift in ground effect by different 
methods, c^ » 2tra; h measured from c/4. 




KEl’KODUClT'lI^V^ 01' mo 

ORIGINAL PAGE IS POOR 




h/c 

Figure 8.- Comparison of computed airfoil lift in ground effect by different 
methods. c^/c = 0.3; 6^ = 20°; a * 2°; = 1.6685; h measured from 





Figure 9.- Comparison of backwash evaluations at control points by different 
methods. y(x) = j/( 1 - x)/x; z/c = 0.05; N = 6. 



.-r 

x/c 


Figure 10.- Comparison of backwash evaluations at vortex points by different 
methods. y(x) * ^(1 - x)/x; z/c * 0.05; N = 6. 
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Figure 11.- Convergence of computed lift coefficient for a rectangular wing 
of AR = 2 oscillating in first bending mode. M = 0.24; k = 0.47. 



Figure 12.- Comparison of computed pressure distributions with experimental 
data for a rectangular wing of AR * 2 oscillating in first bending mode. 
M - 0.24; k » 0.47. 
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(a) Propulsive efficiency. 


.6 



q \ I i i J I I 

0 .1 .2 -3 -4 .5 .6 

k 

(b) Thrust. 

Figure 13.- Theoretical propulsive efficiency and thrust for rectangular 

wings in linear flapping. 
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